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ABSTRACT 

We study the magneto-hydrodynamic tearing instability occurring in a double cur¬ 
rent sheet configuration when a guide field is present. This is investigated by means 
of resistive relativistic magneto-hydrodynamic (RRMHD) simulations. Following the 
dynamics of the double tearing mode (DTM), we are able to compute synthetic syn¬ 
chrotron spectra in the explosive reconnection phase. The pulsar striped wind model 
represents a site where such current sheets are formed, including a guide field. The 
variability of the Crab nebula/pulsar system, seen as flares, can be therefore naturally 
explained by the DTM explosive phase in the striped wind. Our results indicate that 
the Crab GeV flare can be explained by the double tearing mode in the striped wind 
region if the magnetization parameter a is around 10®. 
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1 INTRODUCTION 

In a recent study jBatv et al.ll2013l . also referenced as Paperl 
below), a model has been proposed to address the origin of 
the strong flares observed in X-rays an d gamma-rays from 
pulsars and magnetars environments (IPalmer et al.l l2005l : 
IStriani et al.l l201ll : iTakamoto et all l2014h . These flares are 
short and powerful. Indeed, the observations of Crab pul¬ 
sar nebula show flare duration between a few hours and 
up to severa l days with a rising/falling time scale o f a few 
hours/days (iBuehler et al.l [2012 : IStrimi et ^l2013l l. This 
duration appears to be too large to be associated with the 
rotation of the neutron star (period of 33 ms), and on the 
other hand too short for typical nebula dynamical time scale 
(« 1 year). 

Relativistic magnetic reconnection is usually believed 
to be an efficient mechanism to explain these flares 
lIClausen-Brown fc Lvutikovll20l3 l. possibly in a way sim¬ 
ilar to solar flares where magnetic energy is suddenly re¬ 
leased and converted into other forms of energy. Indeed, the 
highly magnetized wind in the nebula entails a magnetic 
energy reservoir that is large enough to explain the typical 
observed flare energy (« 10®^J). However, achieving a fast 
enough time scale for releasing the magnetic energy is more 
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tricky. For example, the standard view of tearing instabili¬ 
ties developing in a single current sheet leads to a magnetic 
reconnection (with growth of magnetic islands) on a rela¬ 
tively slow time scale that is al so very dependent on the (un¬ 
known) resistive dissip ation (iLvutikov fc Uzdenskvl [2002 ; 
[Komissarov et ahi [200^ . Recently. lUzdenskv et ahl (12011 ) 
and Cerutti et al.l ( 2012h discussed that photons emitted by 
an extremely accelerated pa rticle in an X-l i ne ma y be able to 
explain the Crab flares, and ICerutti et aP ll2014h calculated 
the evolution of a relativistic current sheet in the pulsar wind 
nebula using 3-dimensional PIC simulations. However, they 
did not succeed in reproducing the observed synchrotron 
spectra, in particular the 2011 April flare event. 


In llBatv et al.ll201 j . Paper 1), an explosive mechanism 
is described, that is based on a magnetic reconnection pro¬ 
cess associated to the Double Tearing Mode (DTM). The 
favored site for the emission corresponds to a region situ¬ 
ated in the stripped wind close to the light cylinder radius 
(r « 50r_L). Indeed, the DTM is a Magnetohydrodynamic 
(MHD) instability that is well known to develop in mul¬ 
tiple current sheets, as expected from the magnetic struc¬ 
ture of pulsars magnetosphere with the presence of a cur¬ 
rent sheet wobbling around the equatorial plane (ICoronitil 
Il99(]ll . The explosive character of the mechanism is due to 
the sudden and fast development of a secondary instability 
that is structurally driven by the interaction of the initial 
magnetic islands situated on two successive current layers. 
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It also manifests by a merging of the magnetic islands and 
is characterized by a bulk plasma flow reaching a magnitude 
close to the Alfven speed. 

More recently, the robustness of the mechanism has 
been shown by performing new Relativistic MHD simula¬ 
tions of the DTM of magnetically dominated plasmas, shift¬ 
ing to more realistic high magnetization parameter (ct) and 
magnetic Lundquist number (S which is a me asure of the 
dissip ation) than previously explored in Paper 1 (jPetri et al.l 
|201!tI . also referenced as Paper 2 below). The results indicate 
a weak dependence of the time scale of the explosive phase 
on the magnetic Lundquist number and the magnetization 
parameter, ^0.3 respectively. 

In this paper, with the aim to improve our model, we 
extend the previous results obtained in the two-dimensional 
approximation, by adding the effect of a magnetic guide held 
component perpendicular to the main plane. Such structure 
is indeed expected from a pulsar current sheet structure 
when the rotation and magnetic axis do not coincide. In¬ 
deed, the magnetic held structure in the stripe d wind can b e 
described by an analytical expression given bv IPetril (l2013l) . 
The main magnetic component is toroidal, along the az¬ 
imuthal direction whereas the poloidal part is directed 

along the radius Br- Because of this particular magnetic con- 
hguration, we get B^ oc rn/r and Br oc (rn/r)^ both having 
the same magnitude at the light cylinder. Therefore the ef¬ 
fect of a guide held in the pulsar striped wind should become 
negligible for distances much larger than rn. We re-examine 
the scaling laws dependence with a and S, with a special 
emphasis on the effect of plasmoid-chain, that can develop 
in a transient way in the high S regime. We also investigate 
the energetics of the reconnection process on the basis of 
our approximate magnetized huid model. More precisely, the 
change of the temperature distribution during the explosive 
phase is determined in order to deduce an effective thermal 
particles spectrum, and the ensuing synchrotron emission. 
The results are then compared to the observations. 

The paper is organized as follows. The numerical set-up 
is described in Section (2] The results of the DTM simula¬ 
tions are presented in Section [31 Section [ 3 ] is devoted to the 
computation of the synchrotron spectra deduced from the 
temperature profile during the explosive reconnection event. 
Before ending with the conclusions in Section [51 we discuss 
the implications of our work for the relevant parameters in 
the striped wind, (sheet thickness and other parameters) 
Section [51 


2 NUMERICAL SETUP 

The evolution of the 2-dimensional (2D) double current 
sheets is modeled by using the relativistic resistive magne¬ 
tohydrodynamic (RRMHD) approximation. The RRMHD 
eq uations are solved by using a numerical method developed 
bv lTakamoto fc Inoud (l201ll ) which solves the equations in 
a conservative form using a finite volume method. Magnetic 
field is updated with t he constrained transport algorithm 
llEvans fc Hawle^ll988t) which allows us to treat the diver¬ 
gence free magnetic f ield. In this paper , the numerical set up 
follows the Paper 1 JBatv et al.l 12013 ). with an extra mag¬ 
netic field component Bq added to model the guide field. 

The initial current sheets are described by assuming the 


Table 1. List of the parameters. 


Name 

cr 

S 

Ly 

Ny 

Bg/Bq 

^A,xl ^ 

SAl 

0.2 

200 

10 1 

1024 

0 

0.408 

SA2 

0.4 





0.535 

SA3 

0.6 





0.612 

SA4 

1.2 





0.739 

SA5 

12 





0.961 

SA6 

120 





0.996 

SBl 

0.2 

3200 

101 

2048 

0 

0.408 

SB2 

0.4 


101 

2048 


0.535 

SB3 

0.6 


201 

4096 


0.612 

SB4 

1.2 


201 

4096 


0.739 

SB5 

12 


201 

4096 


0.961 

SB6 

120 


201 

4096 


0.996 

RAl 

0.4 

400 

101 

1024 

0 

0.535 

RA2 


800 


1024 



RA3 


1600 


1024 



RA4 


3200 


2048 



RCl 

12 

400 

10 1 

1024 

0 

0.961 

RC2 


800 


1024 



RC3 


1600 


1024 



RC4 


3200 


2048 



GAl 

12 

1600 

201 

2048 

0.01 

0.961 

GA2 





0.1 

0.956 

GA3 





0.5 

0.859 

GA4 





1 

0.679 

GAS 





1.5 

0.533 

GBl 

0.2 

1600 

201 

2048 

0.5 

0.365 

GB2 

1.2 





0.661 

GB3 

12 





0.859 

GB4 

120 





0.891 

GCl 

12 

200 

201 

2048 


0.859 

GC2 


800 






static relativistic Harris current sheets. The initial magnetic 
field profile is 

B = Bo [1 + tanh[(y - yo)/l] - tanh[(i/ -|- yo)/l]] x -I- Bqz, 

( 1 ) 

and the initial density profile is 

P = Pb +P 0 [l/cosh^[(t/- po)/Z] + 1 / cos}a\{y-yo)/l]] , 

( 2 ) 

where yo = 31 is the half separation of the current sheets, 
and I is the half thickness of the sheets. The gas pressure 
in the sheets is determined by the pressure balance. Note 
also that, for the sake of simplihcation, the guide mag¬ 
netic field component is assumed to be uniform. The ini¬ 
tial temperature T is assumed to be uniform and constant, 
ksT/mP = 1 where kB,m,c are the Boltzmann constant, 
particle mass, and light velocity, respectively. The corre¬ 
sponding sound velocity is Cs ~ 0.516c. The upstream back¬ 
ground density pb is determined by the magnetization pa¬ 
rameter a = B^/Anphc? where h is the specific enthalpy. 
The relativistic ideal equation of state is assumed, in nor¬ 
malized units h = 1 + r/(r — 1)pI pP where P = 4/3. The 
resistivity coefficient rj is assumed to be uniform and con¬ 
stant, and is related to the Lundquist number S through 
S = Ancl/rj (47r is necessary due to the Gauss unit.). Note 
that in this paper, we use the light speed as the character¬ 
istic velocity to define our Lundquist number. 

The simulation domain is bounded by a rectangular do¬ 
main of dimensions [0, L^] x [—Ly,Ly], with a fixed size value 
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Lx = 2QI. As the boundary condition of ^-direction is taken 
to be periodic, the Lx value is precisely chosen in order 
to get a linearly unstable mode with a wavelength close to 
to the fastest one according to MHD stability analysis (see 
Papersl-2 for more detail). The Ly value is taken to be large 
enough to prevent the boundary effects on the sheets evo¬ 
lution, as free boundary conditions are used in this lateral 
^-direction. We divide this simulation domain into homoge¬ 
neous numerical square meshes with size Ly/Ny. 

A small divergence-free magnetic field perturbation is 
also added (at t = 0) as a vector potential 5Az at the cell- 
edges in order to trigger the system evolution: 


5Az = —cos 

ZTT 


^{x-L^/2) 

-L/x 


exp[-(|i/| - yo) ], (3) 


S=200 



where Bi = 3 x 10 “'^Bq is the perturbation amplitude. 

The parameters o,S, Ly, Ny , Be used in the simulations 
are given in Table 1 q The actual numerical value of the 
parameters, such as po,p(,,Bo, follow Paper I. 


3 RESULTS 

As previously reported in Paper 1, the evolution of the sys¬ 
tem exhibits four distinct stages: (1) an initial linear growth 
of two magnetic islands arranged asymmetrically in the two 
current layers and developing on a time scale (2) an 

ensuing saturation on the slow diffusion time scale S^; (3) 
an explosive growth of the islands interacting with a trian¬ 
gular deformation and driving a fast magnetic reconnection 
event characterized by a merging between the two islands; 
(4) a relaxation towards a new stable state. 

This overall behavior does not change in highly magne¬ 
tized plasma regim e with large Lun dquist numbers, as re¬ 
ported in Paper 2 dPetri et al.ll201?^ 'l. An important result 
for the robustness of our model is that, a relatively weak 
dependence of the time scale of the explosive phase on the 
magnetic Lundquist number magneti¬ 

zation parameter (cr°'®) is obtained. 

Interestingly, the explosive stage is however character¬ 
ized by a transient growth of smaller magnetic islands iden¬ 
tified as plasmoid-chain instability. The plasmoids appear 
when a critical local Lundquist number is reached, as it is 
the case for runs using S = 3.2 x 10® (see Fig. 1 in Paper 
2). Thus, in order to quantify the effect of the plasmoid- 
chain, we re-examine in somewhat more details the scaling 
law dependence of the explosive stage. 

3.1 Scaling law dependence and plasmoid chain 
effect 

Following the Papers 1-2, we use the maximum 4-velocity 
in a;-direction, Ux = 714, as a characteristic variable of the 
double tearing mode evolution. 


S=3200 



Figure 1. Maximum 4-velocity yVx of double tearing mode with 
respect to the magnetization parameter a. Top: Low Lundquist 
number cases, S = 200: runs SAl-6. Bottom: Large Lundquist 
number cases, S = 3200: runs SBl-6. 


Firstly, we examine the scaling law of Ux on the mag¬ 
netization parameter for different Lundquist number values. 
The results are shown in Panels of Figure [T] The top and 
bottom panels are for the low and large Lundquist num¬ 
ber cases, respectively. In the low Lundquist number case, 
runs SAl-6, we found that the maximum 4-velocity scales 
as ~ CT®'® and ~ in the small a and large a regimes, 
respectively^. This indicates that the maximum 4-velocity 
in the low-a region is proportional to the Alfven 4-velocity: 
7 aca = However, we note that the obtained maximum 
velocity is not the Alfven velocity but approximately just 
half of it. This is consistent with the outflow velocity of the 
single relativistic tearing instability with small Lundquist 
number, indicating th e deceleration by the large enthalpy 
in the outflow region dTakahashi et al.1 1201 il l. The bottom 
panel of Figure 1 shows the maximum 4-velocity 7I4 in 
the case of S = 3200, runs SBl-6. Differently from the 
low Lundquist number case, the maximum 4-velocity can 
be well reproduced by the Alfven 4-velocity, This is 
because, in the large Lundquist number case, the sheets be- 


^ Note that in this paper, we treated high-Lundquist number 
plasmas. In general, it is very difficult to find the numerical con¬ 
vergence in the presence of plasmoid-chain due to its turbulent 
nature. In this paper, we have checked the convergence of the 
maximum velocity and temperature at the explosive phase of 
DTM, which is given in Appendix A. 


^ Note that the obtained scaling law on tr-parameter is a little 
different from our previous work which indicated ~ This is 

mainly due to including cr = 0.2 which is actually included in a 
new scaling regime we found for the first time in this paper. 


© 2015 RAS, MNRAS OOO.IillTfil 











4 Takamoto, Petri, & Baty 


0=0.4 



t= 11800 I/C, logio(kBT/mc^) 



0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

-0.4 

- 0.6 

- 0.8 

-1 


x/l 


0=12 



100 1000 10000 
S 

Figure 2. Maximum 4-velocity 7 Vx of double tearing mode with 
respect runs RAl-4, SA2. Bottom: High magnetization cases, a = 
12: runs RCl-4, SB5. 


Figure 3. A snapshot of temperature profile at the burst phase 
of run GA4. 


a=12, S=1600 



comes plasmoid-chain as rep orted in Paper 2, which inclu des 
the Petschek type structure llBatvll20ia: iTakamotdli^V 

Next, we examine the scaling law of 7 I 4 on the 
Lundquist number. The numerical results are reported in 
Panels of Figure (2] The top and bottom panels are in the 
small and high a cases, respectively. In the small a case, the 
maximum velocity scales approximately as S'®'® as reported 
in the Paper 1. In the large cr case, the maximum 4-velocity 
scales as ^ g^g small Lundquist number regime. 

However, the scaling becomes ~ S®'® and gradually satu¬ 
rates in the large Lundquist number regime. This is because 
the maximum velocity becomes very close to the Alfven 4- 
velocity in the upstream region, ^/a ~ 3.5, when a — 12. 

3.2 Guide magnetic field effects 

In order to investigate the effect of a magnetic guide field 
component, the direction of the total magnetic field Is varied 
with a fixed magnitude, that is, B = B cos <I>x -f B sin <l?z = 
Bqx + Bqz where B is the constant total magnetic field 
strength, and <I> is the angle between x-axis and the mag¬ 
netic field direction. The corresponding runs are GA1-GA5 
in Table 1. 

Figure [3] is a snapshot of temperature profile of run 
GA4. The basic spatial structure is roughly similar to the 
case without guide field. However, the resulting temperature 
from the burst phase is much smaller than the case without 
guide field (see Figure fTOl) . This is because the compression 


Figure 4. Temporal evolution of the maximum 4-velocity with 
different guide field. 

of the guide field reduces the increase of gas pressure. Note 
also that the sheet region does not have plasmoid-chain, dif¬ 
ferently from the case without guide field as reported in 
Paper 11. This is due to the stabilization o f tearing instabil¬ 
ity by guide field llSomov fc Vernet^ll993ll . which prohibits 
the secondary-tearing instability leading to the plasmoid- 
chain. Figure |4] is the temporal evolution of the maximum 
4-velocity in x-direction. It shows that the initial single tear- 
in g instability is slowed dow n by guide field as is predicted 
by Somov fc Vernete] lll993h in non-relativistic work. It also 
shows the evolution of secondary phase is slowed down and 
the maximum velocity becomes slower as increasing the 
guide field component. 

Figure [ 5 ] plots the guide field dependence of the maxi¬ 
mum 4-velocity. In the figure, the maximum 4-velocity ex- 
hibits a decrease wit h inc reasing guide field as repor ted by 
IZenitani et al.l (l2009l l and lZanotti fc Dumbserj ll201lll . This 
can be interpreted by a decrease of the reconnection mag¬ 
netic field component Bo, which results in the decrease of 
the Alfven velocity, ca = po^l (1 -I- a), along the sheet^. 
This is due to the decrease of the reconnection magnetic 

® Here, = B^/A-wphc^p is the magnetization parameter in 
terms of the reconnection magnetic field component. 
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0=12, S=1600 



Figure 5. The maximum 4-velocity 7 14;: runs GAl-5. 


0=12, S=1600 



Figure 6. The Rutherford time (see text for the definition). The 
horizontal axes are the ratio of the guide field to the reconnecting 
magnetic field component. 


field strength as-. Bq = i3cos<f>. Figure [S] shows the Ruther¬ 
ford time (see Papers 1-2), defined as the time at which 
the explosive phase is triggered. The results show that the 
Rutherford time increases linearly. Note that the slope is 
very steep, and the Rutherford time becomes even twice 
when Bg = 0.5Ro- Note that these properties are very sim¬ 
ilar to the beha vior of the relativ istic Petschek slow shocks 
with guide field llLvubarskvll2005l i . 

Figure [7] are the maximum 4-velocity dependence in¬ 
cluding guide field which is fixed as Bg/Bq — 0.5. The 
top panel shows the magnetization parameter dependence. 
It shows the dependence is completely different from the 
no guide field case in Figure [U and the guide field dras¬ 
tically reduce the maximum 4-velocity, although the maxi¬ 
mum 4-velocity slightly increases with a. This is because the 
guide field hinders the evolution of plasmoid-chain as indi¬ 
cated in Figure [3 and the sheets becomes usual relativistic 
Sweet-Parker sheet whose outflow s aturates an upper limi t 
before reaching the Alfven velocity llTakahashi et al.]|201ll) . 
The bottom panel shows the Lundquist number dependence. 
It shows the maximum 4-velocity increases with Lundquist 
number similarly to no-guide field case shown in Figure [3 
However, the power-law index of the Lundquist number is 
around 0.6 which is an intermediate value indicated in the 



100 1000 
S 

Figure 7. Top: Maximum 4-velocity dependence on cr, run GBl- 
4. Bottom: Maximum 4-velocity dependence on S, GC1,2,GB3. 
In these calculations, the ratio Bq/Bq is fixed. 


Energy Distribution in Reconnected Region 



Bq/Bq 


I I Magnetic Field (B^) I I Thermal 

I I Magnetic Field (Bx.Bz) 1 1 Bulk Kinetic 

Rest Mass 

Figure 8. Energy density in the thermalized region after DTM 
burst with guide field component, corresponding to runs: GA2 
and GA4. 

bottom panel of Figure[2l that is, between 0.4 and 0.8. These 
panels clearly show that the guide field basically weakens the 
energy conversion resulted from DTM, and the burst phase 
becomes less explosive even in highly magnetized plasma 
with high Lundquist number. 

Figure [8] is the energy distribution to each component 
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in the reconnected region of runs GA2 and GA4. It shows 
most magnetic field energy is converted into the thermal 
energy in the both cases. However, the guide field energy 
occupy around 30% of the total energy in GA4. This in¬ 
dicates the guide field magnetic field pressure compensates 
the compression by DTM, and reduces the conversion into 
thermal energ y. This is basically c onsistent with the two- 
fluid result bv IZenitani et ahl (l2009l l. However, the thermal 
energy is still dominant in the case of GA4, and this reflects 
the fact that DTM is not related to slow shocks. 


4 DTM THERMAL ENERGY SPECTRUM 
AND CRAB FLARES 

In our previous works llBatv et al.l l2013l : IPetri et al.l l2015ll , 
we did not discuss in detail the radiation properties associ¬ 
ated to the double tearing mode evolution. In order to obtain 
more quantitative results about possible radiative signatures 
of the reconnection phase, in this section, we compute the 
thermal synchrotron spectrum expected to emanate from 
the relativistic and highly magnetized plasma during its ex¬ 
plosive phase. The procedure is explained in the following 
lines, our aim being to compare our synthetic spectra with 
the variability of the Grab gamma-ray flares. 


4.1 Temperature Profile 

As already discussed in Paper 1,2, the double tearing mode 
shows an explosive reconnection phase during which temper¬ 
ature drastically increases in the reconnected region even in 
the case of moderate value of the magnetization parame¬ 
ter a. In this paper, we explore ultra-relativistic Poynting 
dominated flows assuming a high-Lundquist number for the 
plasmas, and show the resulting explosive reconnection be¬ 
comes more violent, that is, the flow speed becomes faster 
and the temperature becomes higher, as the magnetization 
cr-parameter and the Lundquist number are increased. Since 
the Grab pulsar wind is considered to be a high-a plasma 
with very large Lundquist number, we expect the double 
tearing mode to be sufficiently energetic to explain the flares 
in the Crab pulsar/nebula system. 

Figure [5] shows the profile of temperature just before 
and after the explosive phase of double tearing mode in a 
highly magnetized plasma, a — 120, with high Lundquist 
number, S — 3200. As is indicated in Figure [HI the high 
Lundquist number plasma gives very narrow reconnection 
sheets in which several plasmoids can be observed (see also 
our Paper 2). The plasmoid region becomes very dense and 
also very hot with temperature as high as ksT > amc? be¬ 
cause of the compres sion by nearly A lfvenic reconnection 
flows along the sheet llTakamotdl2013l L After the explosive 
phase, the magnetic field between the two original sheets 
is forced to reconnect, and dissipate its energy into kinetic 
bulk flow and thermal energy. Since the motion inside of 
the reconnected region becomes highly stochastic, the re¬ 
sulting kinetic bulk flow energy rapidly dissipated into the 
thermal energy. In particular, the random motion induces 
strong compression, and this also increase the temperature 
in the plasmoids. 


t = 2250l/c, log.|o (kBl/mc 


10 


-10 


- 2 

- 1.5 

— 1 

- 0.5 

- 0 
-0.5 


20 


x/l 


t = 2300 l/c, log,Q (kg T/mc^) 



2.5 
2 

1.5 
1 

0.5 

0 

0.5 

-1 


x/l 


Figure 9. The temperature profiles logj^Q [/cBT/mc^] during the 
double tearing burst phase of the run SB6. The top and bottom 
panel show just before and after the burst phase, respectively. 


4.2 Synchrotron Energy Spectrum 

Next, we investigate the synchrotron energy spectrum using 
our numerical results. The typical photon energy for the 
synchrotron emission can be written as: 

Z 2 B 2 / A\ 

Esync = 

where B is the magnetic field strength as measured in the 
frame where Csync is detected, Bq — rn}(? jeh is the critical 
magnetic field, m the electron mass, e its electric charge, 7 
the typical Lorentz factor of the electrons and h the reduced 
Planck constant. Assuming the synchrotron radiation is 
mainly emitted by the pair plasma which is in local thermal 
equilibrium due to the MHD approximation, the Lorentz fac¬ 
tor in Equation © can be expressed for an ultra-relativistic 
plasma equation of state as: 7 ~ ZkBT/mc^ = 30. For later 
convenience, we introduced the normalized temperature 0 . 
Note that the temperature T depends on the location in the 
simulation box. In the explosive phase, a strong gradient are 
formed leading to drastic variation in the temperature profile 
from point to point as seen in FigurejO] In the pulsar striped 
wind region, the dominant background magnetic field com¬ 
ponent is toroidal and propagates as an entropy wave thus 
decreasing with radius r according to, Bq = Rl r^/r, where 
rn ~ 1.5 X 10® [m] is the radius of the light cylinder and 
Rl ~ 100[T] the magnetic field strength at rn for the Grab 
pulsar. The synchrotron photon energy in the fluid comoving 
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S=3200 



S=3200 



Figure 10. Top: The energy spectrum with various magnetiza¬ 
tion parameter a. Bottom: maximum temperature with respect to 
background magnetization parameter. In both case, S = 3200 is 
assumed, and the pulsar parameters are r = SOr^, and T\y = 300. 
The runs are corresponded to runs SBl-6, and the data at the 
burst phase reaching maximum 4-velocity is used. 


frame becomes 



where Bq = Bq/Tw is the background magnetic field mea¬ 
sured in the fluid comoving frame, or the simulation frame 
and Fw is the Lorentz factor of the pulsar wind. Finally, the 
Lorentz transformation into the pulsar rest frame (observer 
frame) gives us 

Esync.L = :p; JZ TV = <5w Esync (6) 

Fw (1 — pw cost?) 

where i9 is the angle between the velocity of the reconnecting 
blob and the observer line of sight, /3w is the velocity of the 
wind in the unit of light velocity, and Sw is the “Doppler fac¬ 
tor^’ . If we assume that this blob is pointing towards Earth 
in the most favorable case, dw — 2Fw, the peak energy of 
the photons will be 

£sync,L ~ 0.311[eV] 0^ (l^) 

The energy spectrum can be obtained by integrating Equa¬ 
tion 0 over the numerical spatial domain, f dF'nesync, 


where n is the number density of emitting particles in the 
laboratory frame . The top panel of Figure [10] shows the 
energy spectra with several value for the magnetization pa¬ 
rameter a corresponding to the Runs SBl-6 at the burst 
phase when reaching maximum 4-velocity . Typical values 
for the Crab pulsar are: r = SOrc 0 and Fw = 300 0. This 
panel clearly shows the energy in their body part increases 
with CT-parameter, approximately 100 times increase as a- 
parameter increases by one order of magnitude. This can 
be understood from the bottom panel of Figure 1101 The 
panel is the maximum temperature during the explosive 
phase in terms of the cr-parameter. It shows the maximum 
temperature increases linearly with the cr-parameter. Since 
the photon energy depends on the temperature as T^, as 
indicated in Equation 0, the resulting energy spectrum 
can be estimated as oc cr^. The observed Crab flares 
have their mean energy around 10^ to 10^ [MeVI en ergy re¬ 
gion (iBuehler et al.l[^012l : lBuhler Sz Blandfordl2014l l. so that 
the double tearing mode can explain the Crab flares if the 
Crab pulsar wind has a > 10® which is in agre ement with 
curre nt estimates from theoretical expectations llKirk et al.l 
I 2 OO 9 II about a trans-Alfvenic flow. We also note that the 
energy spectrum in Figure [10] increases with cr-parameter. 
This is because, in high a cases, the n umber of plasm oids 
along the sheets increases as reported in llTakamotcj2013ll ; In 
plasmoid region, the plasma is characterized by a high tem¬ 
perature and a high density, which increases the high energy 
photons and make the energy spectrum harder0 . Note that 
the plasma is far below the radiation reaction limit because 
in the wind frame the energy of the synchrotron photons is 
much less that the limit of 240 MeV. Therefore, the dynam¬ 
ics of the plasma presented in this paper is not significantly 
perturbed by the radiative losses. Our RRMHD code does 
not include any such losses so far. We computed the spectra 
by a post-processing algorithm. 

Finally, the synchrotron cooling time in the wind co¬ 


^ Here, the emission region is assumed to start approximately at 
a fixed ra dius outside the lig ht-cylinder following the prescription 
given by llPetri &: Kirld[2^005h . This is a necessary requirement for 
both fitting the synchrotron spectrum to the observed value (as 
the synchrotron emissivity decreases with distance) and adjusting 
the timescale of the Crab flares taking into account time dilation. 
Although this makes it difficult to discuss timescale of DTM, we 
use the DTM timescale as the observed flare timescale in Equation 
lU since this is one of the characteristic timescales of DTM model 
in the striped wind. The validity of the above assumption will be 
investigated in a future work. 

^ This is twice larger than the value inferred in Paper 2, Pw ^ 
150. This is due to changing of the values of some parameters in 
the dynamical time of the explosive phase, Equation in order 
to explain the cut-off energy of the observed flare. 

® Our numerical results in the case of cr = 120 gave us an ef¬ 
fective index in the high energy region 7 f dx ~ 1-3, which is 
not so far away from the observed energy spectral index of the 
flaring component 7F,obs = 1.27±0.12. Co mparing with the non - 
thermal PIC simulation results obtained bv ICerutti et al] ll20l4) . 
the DTM spectrum is harder in high energy region due to the 
plasmoids. 
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moving frame can be written as: 


^svnr. — - 


4 7(TtC/b 
~2.3 X 10" [s] 

p \ 2 

i W 

300 


0- 


B 


0.05Bo) 




where ctt is the Thomson scattering cross section and Ub 
is the magnetic field energy density in the fluid comoving 
frame. Note that this timescale is estimated in the merged 
sheet where the magnetic field release occurs. In our sim¬ 
ulations, the dynamical time of the explosive phase in the 
simulation frame, fdyn, is about bOl/c. Using I = 27rarw»"L, 
0 it reduces to 



where P = ‘I'ktbIc is the rotation period of the central neu¬ 
tron sta][f|. 

Comparing these 2 timescales, we obtain 

Tsync ^ ^ / £w ^ Q-i f ^ " y 

Tdyn \300y \0.05Boy VSbrny 

( Bb y" f P 

^ Vl00[T]y) V0.05/ V33[ms] 

This indicates the maximum energy by DTM in the Crab 
pulsar wind would be obtained when 0 ~ 10®. Note that we 
obtain a cutoff energy scale around 300 MeV by substituting 
this temperature into Equation © assuming B ~ O.OSBq 
and r = SOrn, which roughly reproduces the observed cutoff 
energy, 350MeV. Figure fTOl also indicates this temperature 
can be obtained when a ^ 10®. 



5 APPLICABILITY TO CRAB PULSAR WIND 

In this paper, we discuss the double tearing mode to ex¬ 
plain the Crab GeV flares. For the simulation setup, we 
assume a separation of two consecutive current sheets to 
be around 3 times of their thickness. Although the sheet 
separation being fairly confidently constrained to be of the 
order of the light cy linder radius using the geometry of the 
striped wind model (lMichellll97ll l. there is still no accepted 
theoretical estimate giving the thickness of one sheet. How¬ 
ever, from an observational point of view, the Fermi Large 
Area Telescope (Fermi/LAT) detected more than hundred 
gamma-ray pulsars showi ng small scale deta ils in their light 
curves in the GeV band llAbdo et al.ll201^ . If we assume 
that the striped wind current sheets are responsible for the 
gamma ray pulses, the width of each pulse should reflect 
the thickness of the sheet itself. Indeed, for sufficiently rela¬ 
tivistic outflows, the spiral structure combined to relativistic 
beaming effects gives rise to pulsed emission. If the current 
sheet has a half thickness of I then the half width of the 

^ o is a ratio between the sheet width and the sheet separation 
(see also Figure [TTJ. 

® Note that if we consider the Lorentz-time dilation this gives a 
flare timescale: Tjyn ~ 2 hours in the pulsar rest frame, which 
is within the observed flare timescale, 8 hours. 


z/rr 



Figure 11. A meridional profile of current sheets in the striped 
wind. The inclination angle is assumed 7r/4, and the light cylinder 
radius rn is used for the unit of each axis. A solid line is z = 
x tan 0 where 6 = tt/A along which the distance between the sheets 
becomes imbalanced. 


corresponding pulse should be IjnfirB where j3 ~ 1. For 
the Crab pulsar, the pulse width is about 10% of the pulsar 
period thus the underlying current sheet thi ckness should 
be around 10% o f the striped wind wavelength llPetri fc Kirlj 
I2005l : lp'etrill2012l l. We consider this supports our MHD treat¬ 
ment of the double current sheets since the indicated sheet 
width is much larger than the kinetic scales, such as the 
electron’s Larmor radius. 

Next, we discuss the conditions required to trigger the 
DTM in the Crab pulsar wind. The first step is the trigger¬ 
ing of the simple tearing mode, which needs kl < 0.6 so that 
I < A/ 10 where A = 2n/k is the wavelength of the pertur- 
batior0. The second step is that the DTM requires kyo < 1 
such that 2nyo < A. In our case j/o = 3Z thus 6ttI < A which is 
approximately comparable to the first condition. According 
to the pulse profile of the Crab, we have I « 27rrL/10 there¬ 
fore simple tearing is triggered for wavelengths A > 27rrL 
which will be easily satisfied since the sheet length in the 
striped wind is much longer than the sheet separation 27rrL. 
The double tearing mode also demands the distance to a 
boundary from sheets should be larger than the sheet sepa¬ 
ration. Although the sheet separation in the striped wind is 
considered to be equal around the equator, in a high latitude 
region, it gradually becomes a combination of long separa¬ 
tion Ai and short separation As which satisfies Aid-As = 27rrL. 
This means the double tearing mode responsible for the GeV 
flare should occur around high-latitude region in the wind as 
indicated in Figure 1111 Since the inclination angle of Crab 

® Although the initial tearing instability needs somewhat long 
time to evolve in hig h-Lundquist number plasma, recent particle - 
in-cell simulations llCerutti et al.l l2015l : Iphilippov et all |2015^ 
show that the current sheets in pulsar winds suffer from strong 
perturbation around light cylinder radius, and we expect such a 
strong perturbation helps the initial tearing instability to grow 
before reaching r = SOrp. 
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Figure 12. A schematic picture of generating guide field in the 
striped wind current sheets via a radial magnetic field compo¬ 
nent in the pulsar magnetosphere. The Green line is the radial 
magnetic field responsible for the guide field in the wind region. 
The blue lines are the poloidal magnetic field responsible for the 
reconnection magnetic field in the wind region. 


0=12, S=800 



t/(Uc) 

Figure Al. Temporal evolutions of the maximum 4-velocity of 
run RC2 using different numerical resolution. 


pulsar is considered to be larg er than 45 degrees (|Rankinl 
I 1993 I : iMoffett fc Hankin3ll999l ). we expect it is always pos¬ 
sible to find a proper latitude in order for the growth of the 
double tearing mode. 

Finally, as is indicated in Figure O the double tearing 
mode becomes less explosive as the guide field increases. We 
consider the guide field in the current sheets will be small 
in the striped wind, and our scenario can still be applica¬ 
ble in the actual Crab pulsar wind. This is because, firstly, 
the radial magnetic held component decreases with increas¬ 
ing the radial coordinate as: Br oc (r/rr)”^, so that it is 
negligible at r ~ 50rL; Secondly, the magnetic held in the 
striped wind is basically generated by the oscillation mo¬ 
tion of the pulsar magnetosphere around the light cylinder 
where the magnetic held is theoretically considered to be 
anti-parallel across the current sheet in the equator. Hence, 
it is likely that the magnetic held in the wind region is also 
nearly anti-parallel across the sheets of striped wind. One 
possibility resulting in a strong guide held component in 
the current sheets is the existence of a monopole-like radial 
magnetic held in the sheet at the light cylinder. As indi¬ 
cated in Figure I12I such a magnetic held will immediately 
change into a poloidal component in the sheet due to the 
oscillation motion, and survive globally in the striped wind 
region. However, the theoretical studies of pulsar magne¬ 
tospheres do not support this monopole-like magnetic held 
structure, and this will be unlikely to occur. For these rea¬ 
sons, we expect the guide held to be very small, and does 
not reduce the energy release. 


6 CONCLUSIONS 

The striped wind structure represents a natural magnetic 
held conhguration arising from the rotation of a magnetized 
neutron star surrounded by a relativistic pair plasma. It gen¬ 
erates a current sheet wobbling around the equatorial plane. 
From simple analytical models, it is known that the mag¬ 
netic held lines become mainly toroidal far from the neutron 
star with a poloidal component decreasing much faster than 
the toroidal component. Therefore, locally, the striped wind 
can be depicted as a current sheet with a guide held for 


which the strength depends on the distance from the light- 
cylinder. 

Our new RRMHD simulations of the DTM including 
guide held effect conhrm our previous works. DTM is a good 
candidate to explain the variability of the gamma-ray emis¬ 
sion of the Crab flares. Indeed the explosive phase exhibits a 
magnetic reconnection event releasing magnetic energy into 
heat and therefore also into synchrotron radiation. The com¬ 
parison between the synchrotron cooling time and the dy¬ 
namical timescale shows that the plasma regime is far be¬ 
low the radiation reaction limit. The spectrum peaks at a 
energy Sp simply related to the magnetization parameter a 
by Ep oc a^. This allows us to constrain the magnetization 
of the Crab to lie around cr « 10®. This is in addition to 
the previous constraints on the location of the origin of the 
flares estimated to be around r « 50 rr where rr = c/Q 
Is the light-cylinder radius, c the speed of light and D the 
rotation speed of the pulsar and on the maximum L orentz 
factor of the striped wind, F < 300 (IPetri et al.ll201^ . 
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APPENDIX A: CONVERGENCE OF THE 
CALCULATIONS 

In this paper, we dealt with high-Lundquist number plas¬ 
mas. In general, high resolution is necessary to simulate 
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high-Lundquist number plasma to reduce numerical dissipa- This paper has been typeset from a T 13 X/ file prepared 

tion. Figure [XT] shows temporal evolutions of the maximum by the author, 

velocity of run RC2 with different resolutions. It shows the 
resolution approximately Nx > S/2 is sufficient for repro¬ 
ducing the maximum 4-velocity. Note that our code does 
not include viscosity and thermal conductivity but only re¬ 
sistivity, so that it is in general impossible to obtain perfect 
numerical convergence. 
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